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Abstract. Identified particle observables from viscous hydrodynamics are 
sensitive to the fluid-to-particle conversion. Instead of the commonly assumed 
"democratic" Grad ansatz for phase space corrections 5f, wc utilize corrections 
calculated from linearized covariant transport theory. Estimates based on a tt — p 
system with binary collisions indicate that protons are much closer to equilibrium 
than pions, significantly affecting the dissipative reduction of differential elliptic 
Qs 1 flow in Au-|-Au at RHIC. In addition, we test linear response against fully 

*>/| ' nonlinear transport for a two-component massless system in a Bjorken scenario. 

Strikingly, we find that, while linear response accounts well for the dynamical 
sharing of shear stress, the momentum dependence of phase space corrections is 
best described by Grad's quadratic ansatz, and not the linear response solution. 
-)— » , 

^ ' PACS numbers: 24.10.Lx, 24.10.Nz, 25.75.Ld 
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1. Introduction 

^ I An inevitable ingredient in contrasting hydrodynamic calculations of heavy-ion 

^-^ ■ collisions to experimental data is the conversion of the fluid to particles (hadrons). 

For ideal fluids this is straightforward, at least within the Cooper- Frye framework [I]. 
(y-s ■ For viscous fluids, however, dissipative corrections such as shear stress, distort phase 

space distributions from local equilibrium, / = feq + 6f, and the very same viscous 
^— V ■ hydrodynamic fields can be described by an infinite class of phase space corrections. 

Viscous hydrodynamic calculations[2l[3l|4] commonly ignore this ambiguity and, 
in the presence of shear stress, pick corrections for each hadron species in the 
"democratic" form ([1]). This ad- hoc choice does not reflect that species that interact 
more frequently should be closer to equilibrium. 

We use here covariant transport theory to compute dissipative phasespace 
S . corrections in a hadron gas. In the linear response limit, quite akin to [S], we 

obtain first results for a pion-proton system with binary 2 — )■ 2 rates, illustrating how 
dynamical sharing of shear stress between species affects identified particle differential 
elliptic fiow V2{pt)- We also test the validity of linear response against fully nonlinear 
2 — > 2 transport for a two-component mixture undergoing a 0+lD Bjorken expansion. 

2. Dissipative phase space corrections and covariant transport 

For zero bulk viscosity, the only dissipative field in viscous hydrodynamics is shear 
stress tt'"^ (x) and therefore dissipative corrections for species i are of the form Sfi = 
Ciipiau" /T)'K'^'^Pif^Pi^fi^eq, whcrc u IS thc flow velocity (Boltzmann statistics used). 
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Ci are arbitrary functions with the only constraint that partial shear stresses sum to 
the total Y^ 7r[' 



n^" , naturally satisfied by the "democratic" Grad prescription 

1 , ttA"'- 
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2(e + p)r2 ■'* 2r2(e+p) 

On the other hand, from covariant transport one expects 5fi to depend on the 
mean free paths, or equivalently the inverse Knudsen numbers Ki ^ 'Ynj{avrei)ij /t. 
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We consider on-shell[3] 2 — ?► 2 Boltzmann transport for each species. 

Linear response to shear is the same expansion around local equilibrium as the 
standard calculation of shear viscosity [71 18], leading to an integral equation for {Sfi} 
that is solved variationally, resulting in corrections proportional to spatial gradients 
in the comoving frame (as in Navier-Stokes hydrodynamics). Because nonequilihrium 
terms are only kept in the collision kernel, this ignores shear stress relaxation captured 
by second-order hydrodynamics P [TOl [H [31 [S]. 

Alternatively, we solve the fully nonlinear transport for a two-component mixture 
that evolves near equilibrium. This way we study 5f in a less constrained dynamical 
setting, which also provides useful cross-checks for the linear response results. 
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Figure 1. Differential elliptic flow V2{pt) of protons (symbols) and pions (lines 
witli no symbols) for various dissipative corrections Sf and freezeout temperatures 
Tfo = 165 MeV (left) and 140 MeV (right). Dotted: ideal ffuid (<5/ = 0), dashed: 
viscous fluid with "democratic" Sf, solid: viscous fluid with Sf determined 
dynamically for a tt — p mixture. In both viscous cases, we take the Grad ansatz 
Sf <x p^ and set the same partial shear stress for pions (i.e., the solid and dashed 
lines without symbols coincide). The tt — p splitting pattern is sensitive to the tt 
and p equilibration rates, especially at higher Tf^. 



3. Highlighted results 

Figure [T] shows how dynamical shear stress sharing influences differential elliptic flow 
V2{pt) at freezeout in Au-f Au at RHIC. We focus here on a tt — p system, i.e., ignore 
interactions of protons and pions with other hadrons, and also ignore resonance decays. 
The "dynamical" results are from linear response with isotropic, energy-independent 
effective scattering cross sections o-jrTr = 30 mb, (T,rp = 20 mb, app = 50 mb. At 
T ~ 120 — 165 MeV, these approximate quite well a more realistic calculation [11] of 
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mean collision frequencies in a tt — p system based on measured phase shifts. For 
simplicity, quadratic momentum dependence was assumed for both species, however, 
C^ and Cp were optimized variationally, i.e., C^ j^ Cp. The "democratic" prescription 
corresponds to quadratic momentum dependence with C,r = Cp. Results for ideal 
hydrodynamics {Sf = 0) are also shown. Ideal hydrodynamic fields for Au+Au at 
^/s = 2Q0A GeV, 6 = 7 fm, were obtained with our patched version 0.2p2 of the 
2+lD boost-invariant AZHYDRO code[12]. We used the recent s95-pl equation of 
state parameterization[13] matching lattice QCD results to a hadron resonance gas. 
Shear stress at freezeout was estimated as in[T3], based on gradient corrections to ideal 
hydro tt'^'" = 7/[V^m'' + V''u^-(2/3)A^'^(9aW")] with A'^'^ = ^^'^-u^m", V = A'^'u^. 
But unlike [14], we use real hydro instead of the parameterizations ( "Blast- wave" ) . 

The left plot in Figure [T] shows the tt — p splitting of V2{pt) for freezeout 
at T — 165 MeV, the typical switching temperature between hydrodynamics and 
transport in hydro-|-transport calculations [TS]. At moderate pT ^ 1 — 2 GeV, we 
find a significant viscous reduction of elliptic flow for both species even for a small 
shear viscosity to entropy density ratio rj/s — 0.1. However, compared with the 
dynamical approach in which protons have smaller viscous correction than pions, the 
"democratic" prescription oversuppresses proton elliptic flow by 20 — 30%. In the 
dynamical approach, pion-proton crossing occurs at much lower pT ^ 1 GeV. At 
lower T = 140 MeV (right plot), more applicable to pure hydrodynamic description 
of RHIC data, dissipative corrections are smaller, but the relative difference between 
the "dynamical" and "democratic" approaches is in fact larger. 

Figure [2] shows dissipative corrections from fully nonlinear 2-^2 transport for a 
two-component massless gas A -\- B undergoing 0-l-lD Bjorken expansion (transverse 
translational symmetry assumed) . /a and fs only depend on normalized proper time 
t/tq^ Pt, and the difference between coordinate and momentum rapidity. We set 
isotropic oaa '■ o'ab ■ cfbb =4:2:1, densities ua = ns, such that (initial) inverse 
Knudsen numbers are Ka{to) — 2, Kb{to) = 1, i.e., species A scatters twice as often 
as B, and start from local equilibrium. To keep r]/s « const, we scale[6] cross sections 
with time as a a t^/^. We test how well various ansdtze for Sf can reconstruct, solely 
from the temperature and partial shear stresses, the spectra from the transport. 

With the dynamically determined shear stresses from the transport, Grad's 
quadratic ansatz Sf oc p^ (left plot) is accurate to within about a percent (!) for a 
wide range of times (up to IOtq), and momenta (up to 5 — 6 x T). However, with the 
momentum dependence 6f ex p^ calculated from linear response (right plot), the 
agreement is much poorer. The likely reason is that shear stress relaxation competes 
with the rapid expansion of the system, and therefore, relaxation to Navier-Stokes 
should not be assumed. Note, the power 1.47 we obtained for isotropic 2—^2 processes 
is quite similar to the one in [F for forward-peaked 2-^-2 and radiative 1 <->■ 2. 

4. Conclusions 

Reliable extraction of shear viscosity from heavy-ion data must involve proper 
conversion of viscous fluid to hadrons. Though several aspects of this work need to be 
refined in the future (e.g., with realistic cross sections, real viscous hydro, quantum 
statistics, more species), our results do indicate that the widely used conversion 
prescription ("democratic" Grad approach) is inconsistent with scattering rates in 
a hadron gas. This brings into question the accuracy of viscous hydrodynamic 
calculations of identified particle observables to date. 
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Figure 2. Ratio of spectra from covariant transport to reconstructed spectra 
using only the partial shear stresses from the transport and an ansatz for 5/, 
for a two-component massless system in a 0+lD Bjorken scenario. The ratio 
is plotted at various proper times r = const X tq versus normalized transverse 
momentum Pt/T{t), for species A, which is set to scatter twice as often as B (the 
ratio of inverse Knudsen numbers is Kj^/ Kg = 2). We find that Grad's quadratic 
ansatz for the momentum dependence (left plot) is much more accurate than the 
momentum dependence computed from linear response (right plot). 
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